Possible triplet superconductivity in MOSFETs 
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A theory that predicts a spin-triplet, even-parity superconducting ground state in two-dimensional 
electron systems is re-analyzed in the light of recent experiments showing a possible insulator-to- 
conductor transition in such systems. It is shown that the observations are consistent with such an 
exotic superconductivity mechanism, and predictions are made for experiments that would further 
corroborate or refute this proposal. 



Since 1979 it had been generally believed that two- 
dimensional (2-d) electron systems cannot undergo a true 
metal-insulator transition (MIT), but rather are always 
insulating at zero temperature (T = 0), even for arbi- 
trarily weak disorder. Inititially, this was the-.conclusion 
only for modglp, of noninteracting electronsjd but later 
developmentsB □ strongly suggested that it remains valid 
in the presence of electron-electron interactions. Recent 
experiments on Si^MOSFETs have challenged this con- 
ventional wisdom.ErS The samples used in these experi- 
ments are in an unprecedented parameter regime, as they 
achieve high electron mobilities at low electron densities. 
The apparent quantum phase transition from an insulat- 
ing phase to a conducting one that is observed in these 
samples occurs at very low densities, which leads to a 
very strong effective electron-electron interaction. The 
disorder is also very strong, as is indicated by a 'sepa- 
ratrix' value of the sheet resistance that is about three 
times the Mott number. A MIT has also been observed 
in p-type SiGe quantum wells with paramjeter values that 
are very close to those in the MOSFETs.B 

These observations can in principle be theoretically ex- 
plained in a number of ways. For example, it is known 
that there is a MIT in 2-d systems with spin-orbit scat- 
tering and either no interactioa, or a short-ranged in- 
teraction between the electrons.Q While this universality 
class has beea invoked in the context of the experiments 
in question,E3 it is unclear why it should apply to the 
systems under consideration, where the electrons inter- 
act by means of a strong, long-ranged Coulomb inter- 
action. Other possible theoretical explanations include 
superconductor-to-insulator tMinsitions (SITs), which are 
well known to occur in c? = 2.eI Indeed, it has been noted 
in Refs. ^ and [t] that the observed transition has many 
features in common with SITs. It has also been pointed 
out in Refs. pT|and[T^ that the conducting phase is likely 
to be a superconductor, or a perfect conductor, respec- 
tively. 

Some time ago the present authors have discussed 
a mechanism that can lead to spin-triplet, even-parity 
superconductivity in disordered, interacting electron 



systems. El The predicted effect is strongest in d = 2, and 
MOSFETs were explicitly mentioned in Ref. 13 as po- 
tential realizations of this phenomenon. This motivates 
us to analyze this mechanism semi-quantitatively for pa- 
rameter values that are appropriate for the new MOS- 
FET samples. As we will show below, the observations 
are compatible with an explanation in terms of super- 
conducting fluctuations that preceed the predicted ex- 
otic superconducting state, and we are also able to make 
predictions for observables other than the resistivity. 

The physical mechanism that underlies the exotic su- 
perconductivity discussed in Ref. |l3| is based on the 
slow decay of charge and spin polarization clouds that 
is caused by long-time tail effects in disordered systems. 
The details have been explained in Ref. |l^, and need not 
be repeated here. We just recall that there are two physi- 
cally distinct mechanisms, based on charge and spin fluc- 
tuations, respectively. Characteristic features of the en- 
suing even-parity, spin-triplet superconductivity are, (1) 
a very low Tc for parameter values that are characteris- 
tic of conventional 2-d electron systems, (2) a very rapid 
variation of as a function of parameter values such as 
the bare interaction amplitudes and the disorder, (3) a 
gapless tunneling density of states, and (4) a logarithmic 
dependence of the specific heat coefficient on the temper- 
ature. All of these features were derived and discussed 
in Ref. |l3|. In addition, one finds (5) an ordinary Meiss- 
ner effect, like in a gapless BCS superconductor, as we 
will now proceed to show. To this end, we calculate the 
transverse current-current susceptibility, using the field- 
theoretic formulation of Ref. |l^. A standard calculation 
at the same level as our previous saddle-point theory for 
the order parameter yields a transverse current suscepti- 
bility 
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Here Q and J- are the normal and anomalous Green func- 
tions, respectively. 
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with LOn a fermionic Matsubara frequency, = p^/2to — 
Ei? with ep the Fermi energy, and A„ the frequency de- 
pendent gap function. The gap function on the imaginary 
axis is purely real. An = — A_„ = A(ia;„) = (5 (a;,,), with 
i5(ix') a real function.EJ At zero temperature, Eq. ( [la] ) can 
then be expressed as a frequency integral, 
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with n the electron number density. We see that the 
correction to the normal- metal result, n/m, is negative 
definite.Ej Consequently, we have XT(k — + 0) < n/m, 
which guarantees ideal diamagnetism. 

Si MOSFETs were discussed in Ref. ^ as possible re- 
alizations of the proposed mechanism, but for typical pa- 
rameter values at best ultralow values for Tc, in the fiK 
region, were to be expected. For the samples investigated 
in Refs. |^-^, the situation is very different. Let us esti- 
mate the mean-field Tc for these samples, assuming the 
charge fluctuation mechanism described in Ref. g_3| (we 
will come back to the spin fluctuation mechanism later). 
The mean-field value of Tc is given by 



Tc = Tf ^ cxp 
Rn 
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Here n — k{1 + Fq), with k the Thomas-Fermi screening 
wavenumber, and Fq a Landau parameter. Tf is the 
Fermi temperature, and Rn is the resistance per square 
measured in units of h/e'^ = 4108ri. With an electron 
density n = 10^^ cm"^ and a valley degeneracy gv — 2, 
we have a Fermi wavenumber kp = 5.6 x 10^ cm~^. With 
an effective mass m* = 0.2 me, this corresponds to a 
Fermi temperature Tf = 6.9 K. Assuming a Landau 
parameter Fq — —0.9, the screening wavenumber is k = 
1.5 X 10^ cm^^. At a sheet resistance equal to i? — 
0.8 X 27r?i/e2, the Tc formula, Eq. (|), yields Tc w 1.7 x 
lO-^Ti. « 10 mK. 

In the experimental range of Refs. one therefore 
does not expect to see true superconductivity resulting 
from the charge fluctuation mechanism. If the observed 
transport anomalies are to be explained in terms of it, 
then they must be caused by superconducting fluctua- 
tions. Let us therefore estimate the width of the temper- 
ature region within which one expects substantial fluctu- 
ations. To this end we recall that the kernel in the gap 
or Tc equation in Ref. |l3| takes the form of a function F, 
whose argument depends on T and Ra- For asymptot- 
ically small arguments, F approaches a logarithm, and 
we have 



F{Ru{kF/RfT/TF) - \n{Ru{kF/RfT/TF) = f{T). 

(3) 

If we consider T the flow parameter in a RG description 
of the system, then Ru becomes T-dependent, and for 
fixed bare parameters F can be considered a function 
/ of T, as is indicated in Eq. (^). If we replace T by 
Tc, and use the above parameter values, then we obtain 
f{Tc) ~ —10. As a criterion for noticeable fluctuations, 
we require the temperature to be lower than some T* de- 
fined by f{T*) = — 5. As the temperature increases from 
Tc, the scale dependent quantity Ra, increases, but in 
d = 2 it does so only logarithmically.H The most impor- 
tant T-dependence is therefore the explicit one in Eq. (|^). 
Keeping only this leading T-dependence, the estimate be- 
comes structurally the same as in BCS theory. Using the 
same parameters as for our estimate of Tc above, we find 
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We also have obtained an estimate of the same order by, 
using the RG flow equation for the disorder parameter 
which determines the T-dependence of the resistivity. 

With the above values for Tc and T* , in a 2-d sys- 
tem with parameter values appropriate for the new MOS- 
FETs, Eq. (^ predicts that observable superconducting 
fluctuations are to be expected up to temperatures of 
about IK, and the fluctuations should be substantial at 
the lowest temperatures reached in the experiments. This 
could well explain the observations. We stress that these 
considerations amount at best to an order-of-magnitude 
estimate, while reliable calculations of the superconduct- 
ing Tc are notoriously difficult for any mechanism. Nev- 
ertheless, the results show that our mechanism is a viable 
candidate for explaining the observations. It is also en- 
couraging that our Tc is an extremely rapidly varying 
function of the parameter values. This may explain why 
similar effects were not observed in other 2-d systems. 

A magnetic field weakens our superconductivity mech- 
anism, since it quenches the strong renormalization of 
interaction parameters that leads to the pairing. In 
Ref. |l^ we estimated the upper critical field to be 
Hc2 ~ ific/e)kpTc/TF. At the fluctuation tempera- 
ture T* > Tc, one would then expect a field of strength 
H* = {hc/e)kpT* /Tf to have a qualitative effect on 
the transport properties. At T* = 1 K, this yields 
H* — 9 kOe. This is in very good agreement with the 
observed S|UBpression of the conducting phase by a mag- 
netic fieldJl3 

While the above estimates show that our proposed 
mechanism is definitely compatible with the existing ex- 
periments, other observations are needed to convincingly 
corroborate or refute our proposal. Of particular interest 
are the properties of the superconducting phase, if any, 
that we predict to exist at sufficiently low temperatures. 
Apart from the prediction of a Meissner effect, we have 
calculated three observables, viz. the gap function and 
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FIG. 1. Real part (solid line) and imaginary part (dashed 
line) of the gap function A at T = as a function of the 
frequency u, the tunneling density of states A'^ at T = as a 
function of the frequency, and the specific heat coefBcient 7 
as a function of the temperature T as predicted by the theory 
for the parameter values discussed in the text. A'^ and 7 are 
normalized by the free electron density of states per spin, Nf, 
and A, oj, and T are measured in Kelvin. 



the tunneling density of states, both at T = as a func- 
tion of the frequency, and the specific heat coefficient, 
7 — Cv/T, as a function of T, in the superconduct- 
ing phase. Such calculations for any purely electronic 
superconductivity mechanism pose an intrinsic problem, 
namely that the feedback of the superconducting statp.Qn. 
the coupling mechanism must be taken into account £II'E£I 
This problem has so far not been solved. In the case 
of the present mechanism it is likely to just lead to a 
smaller effective coupling constant A = -Rq/tt^, i.e. a 
smaller Ra , then the one that enters the Tc equation, 
except for a qualitative effect in the asymptotically low 
frequency regime £3 For our present illustrational pur- 
poses we choose an effective coupling constant A = 0.4/7r, 
which is a factor of 4 smaller than what corresponds to 
the Ru used for the Tc estimate above. The results are 
shown in Fig. |l|. Notice that the purely electronic pairing 
mechanism leads to a pronounced pseudo-gap structure 
on a frequency scale that is, for our value of the cou- 
pling constant, some 5,000 times larger than Tc- The Tc 
scale is reflected in the region of validity for the asvaap- 
totic low- frequency and low-temperature behavior .lij It 
turns out that the overall shapes of the curves do not 



strongly depend on the effective A chosen, only the en- 
ergy scales, do. The two relevant energy or frequency 
scales areB 17* = (2 x IQ-"^ /Xj^K/kpfTp, which sets the 
overall frequency scale, and = exp(— IZ-s/A), which 
separates the low and high-frequency regimes (Jl* 2 are 
the frequency scales wjf 2 of Ref. multiplied by tt/IOO). 
In conventional, phonon-induced, superconductivity the 
role of Vl\ is played by the Debye frequency, and that of 
by the value of the gap function at zero frequency. 
For the parameter values chosen, the numerical values of 
the two frequency scales are Q,\ — 778 K, and — 47 K. 

Figure |l| shows our qualitative predictions for the su- 
perconducting state. For experiments above Tc, the most 
promising quantity to measure would probably be the 
tunneling density of states, which should reveal the de- 
velopment of the characteristic structure that is shown 
in Fig. |l|. Although the full features would be present 
only in the superconducting state, one would expect some 
structure in the DOS in the region where the transport 
properties are strongly anomalous. This structure should 
become stronger with decreasing temperature, but not 
develop into a real gap, in contrast to a BCS supercon- 
ductor, but rather into the pseudogap shown in Fig. |l|. 

We now briefly discuss the spin fluctuation mecha- 
nism, which in principle is capable of producing much 
higher values of Tc than the charge fhojptuation mecha- 
nism. However, for technical reasonscJ it is harder to 
make semi-quantitative predictions in this case. The Tc- 
formula derived in Ref. O is 



Tc = To exp 



32 



(5) 



Here 7^* is a bare (i.e., high-temperature) spin-triplet in- 
teraction constant, and Tq is a microscopic temperature 
scale that is on the order of the Fermi temperature. If we 
assume 7^ = 1 , which corresponds to a high-temperature 
spin susceptibility that is twice that of free fermions with 
mass m*, and the disorder as above, then we obtain a 
value for Tc that is a substantial fraction of Tp. This is 
certainly an overestimate, for reasons discussed in Ref. 
p^ . However, the point is that the spin fluctuation mech- 
anism allows for much higher Tc-values than the charge 
fluctuation mechanism, provided that there are strong 
spin fluctuations in the system, i.e., a large magnetic 
susceptibility in the normal metal. If the spin suscep- 
tibility of the 2-d electrons in MOSFETs could be mea- 
sured, then it would be very interesting to look for a 
correlation between an enhanced spin susceptibility and 
the occurrence of the metallic phase. The qualitative fea- 
tures of the other obervables discussed above are virtually 
unchanged for the spin fluctuation mechanism. Explicit 
calculationgli3 of the gap function, the density of states, 
and the specific heat, yield results very similar to those 
shown in Fig. |l|, but with the characteristic features oc- 
curring on a smaller frequency or temperature scale. 

We conclude by means of two additional remarks. 
First, there has been a recent report about a similar 
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conductor-to- insulator transition in Si/SiGe heterostruc- 
tures with a much smaller valueJor the disorder than in 
the MOSFETs discussed above.li^ However, the interpre- 
tation of these data has been called into questionJ^J and 
more experimental work seems to be needed in order to 
ascertain if the phenomenon indeed persists to such low 
disorder. It is worthwhile to point out that long-time 
tails exist even in clean Fermi liquids, and a mechanism 
for triplet superconductivity similar to the one discussed 
above is conceivable. The effect would be weaker, though, 
and it would probably not result in a superconducting 
ground state even in d = 2. 

Second, there are other possible interpretations of the 
experiments, which are not related to superconductiv- 
ity. For instance, it was pointed out in Ref. |l| that the 
RG flow equations for disordered interacting electrons in 
d — 2 contain a nonperturbative critical fixed point that 
corresponds to a MIT. Since the fixed point is separated 
from the weak disorder regime by a region of runaway 
flow, its physical signifioajice was considered doubtful. 
However, Castellani et al.Ej have recently speculated that 
this flxed point might describe a MIT in d = 2. In this 
context we point out that one way for the flxed point 
to attain a physical meaning would be to make physical 
sense out of the runaway flow, e.g. |-ip terms of a non- 
Fermi liquid state at small disorder.cj Indeed, in d = 3, 
these runa"\piy flows have been related to the ferromag- 
netic stateo The authors of Ref. ^ have also argued 
that everLJii the absence of a true MIT, the known flow 
equationaJ'tl for disordered interacting electrons in d = 2 
allow for metallic behavior over a sufficiently large tem- 
perature region to explain the experimental observations. 
This would require, however, a large spin susceptibility 
in the metallic phase, large enough to trigger supercon- 
ductivity via the spin-triplet mechanism according to our 
estimates. Thus a combination of some of these sugges- 
tions is conceivable: The pre-asymptotic analysis of Ref. 
^ could explain the existing experimental data, but at 
lower temperatures a transition to the superconducting 
state discussed in this paper could occur. In either case, 
experiments looking for an enhancement of the magnetic 
susceptibility would be very interesting. 
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